arXivil509.05627v2 [quant-ph] 4 Jan 2016 


Transport of ultracold atoms between concentric 
traps via spatial adiabatic passage 

J. Polo^, A. Benseny^, Th. Busch^, V. Ahufinger^, and 
J. Mompart^ 

^ Departament de Fisica, Universitat Autonoma de Barcelona, E-08193 Bellaterra, 
Spain 

^ Quantum Systems Unit, Okinawa Institute of Science and Technology Graduate 
University, Okinawa, 904-0495, Japan 

Abstract. Spatial adiabatic passage processes for ultracold atoms trapped in tunnel- 
coupled cylindrically symmetric concentric potentials are investigated. Specifically, we 
discuss the matter-wave analogue of the rapid adiabatic passage (RAP) technique for a 
high fidelity and robust loading of a single atom into a harmonic ring potential from a 
harmonic trap, and for its transport between two concentric rings. We also consider a 
system of three concentric rings and investigate the transport of a single atom between 
the innermost and the outermost rings making use of the matter-wave analogue of the 
stimulated Raman adiabatic passage (STIRAP) technique. We describe the RAP-like 
and STIRAP-like dynamics by means of a two- and a three-state models, respectively, 
obtaining good agreement with the numerical simulations of the corresponding two- 
dimensional Schrodinger equation. 
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1. Introduction 

Toroidal trapping potentials [1, 2, 3, 4, 5, 6] for ultracold atoms are essential in many 
quantum metrology implementations, e.g., for matter-wave Sagnac interferometry [7], 
and as elementary building blocks in the emerging field of atomtronics [8]. In fact, 
ring potentials constitute the simplest nontrivial closed-loop circuits and offer unique 
features to investigate, for instance, superfluidity and persistent currents [9]. Angular 
momentum can be transferred to the atoms in such systems using either light [10] or by 
rotating a tunable weak link [11]. The latter has been one of the key elements that has 
opened the possibility to study matter-wave analogues of the superconducting quantum 
interference devices (SQUIDs) [12], Recently, hysteresis, which is a basic ingredient in 
electronics, has been observed experimentally between quantized circulation states of a 
Bose-Einstein condensate (BEC) trapped in a ring [13]. 

BECs in two tunnel-coupled concentric ring potentials have also been investigated, 
paying special attention to the generation of topological phases under rotation [14]. 
Ground-state phases, produced by the interplay between two-body collisions and spin- 
orbit coupling [15] or between the intra-atomic and inter- atomic interactions in a 
two-component BEC [16], have also been analyzed. In this context, the development 
of techniques allowing for a high hdelity and robust transport of ultracold atoms 
between ring potentials is a topical issue. Spatial adiabatic passage (SAP) processes 
[17] have been proposed for the transport of single atoms [18, 19, 20, 21, 22], electrons 
[23], and Bose-Einstein condensates [24] between the outermost traps of in-line triple¬ 
well potentials. Moreover, two-dimensional (2D) SAP processes for three tunnel- 
coupled harmonic wells in a triangular conhguration have recently been discussed 
for matter-wave interferometry [25], and for the generation of atomic states carrying 
orbital angular momentum [26]. It is also worth to highlight that SAP has already 
been experimentally demonstrated for light beams propagating in three evanescently 
coupled optical waveguides [27], and recently discussed for the manipulation of sound 
propagation in sonic crystals [28]. 

In this work, we investigate the use of 2D SAP processes for the transport 
of ultracold atoms between the localized ground states of cylindrically symmetric 
concentric potentials. Specihcally, we focus on the matter-wave analogues of the well- 
known quantum-optical rapid adiabatic passage (RAP) [29] and stimulated Raman 
adiabatic passage (STIRAP) [30] techniques. We consider the matter-wave analogue 
of the RAP technique for the loading of a single atom from a harmonic trap into a ring 
potential, see Fig. 1(a), or for transporting it between two concentric ring potentials, see 
Fig. 1(b). For these RAP-like processes, we vary the trapping frequency and the radius 
of the outer ring potential in order to allow for adjustments of the tunneling rate and 
the energy bias in such a way that the atom adiabatically follows a particular energy 
eigenstate of the system. Moreover, by means of a STIRAP-like technique, we investigate 
the transport of a single atom between the innermost and outermost traps of a triple¬ 
ring potential, see Fig. 1(c). In this case, we assume equal radial trapping frequencies 
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for the three harmonic rings and the tunneling-based transport is achieved by modifying 
the rings separations. Both SAP processes studied here rely on the ability to vary the 
frequencies and radii of the trapping potentials. For optically generated potentials, the 
frequency can be modihed by varying the laser intensity that produces the potential [2] , 
and the radius of the ring can be changed dynamically using time-averaged adiabatic 
techniques [3, 4]. These matter-wave analogues of RAP and STIRAP are described with 
a two- and a three-state model, respectively, and the obtained results are found to be 
in good agreement with the corresponding numerical simulations of the 2D Schrddinger 
equation. Finally, we also investigate the possibility to use SAP processes between three 
concentric ring potentials to transport orbital angular momentum (0AM) states. Such 
transport can play a very important role in future atomtronic devices, where 0AM states 
are one of the building blocks to study SQUIDs [12]. However, angular momentum states 
also introduce extra degrees of freedom which consequently will affect the adiabaticity 
condition of the 2D SAP processes. It is therefore necessary to check if the adiabatic 
techniques are still usable in ring type geometries. For instance, in ID systems, the 
energy biases are only affected by the difference in frequencies of the traps. In 2D SAP 
systems, however, the energies of the rings also depend on their radii and thus, coupling 
and detuning cannot be treated as independent parameters. 

The article is structured as follows. The physical system under investigation 
is introduced in Section 2. Section 3 is devoted to the transport of a single atom 
between two concentric potentials by means of the matter-wave RAP technique, whereas 
Section 4 focuses on the transport of an atom between the innermost and outermost 
traps of a triple-ring potential through the matter-wave STIRAP technique. Section 5 
presents the conclusions. 


2. Physical system 


We consider a single atom trapped in 2D cylindrically symmetric external potentials 
formed by different concentric combinations of harmonic and harmonic ring potentials, 
see Fig. 1. Using cylindrical coordinates (r, 0), the harmonic potential with trapping 
frequency Uh for an atom of mass m reads 


Vh{r) = 


and its ground state wavefunction is 

f 


TflUJh ^2 

- P 2h ' 

TTh 


( 1 ) 

( 2 ) 


The harmonic ring potentials of trapping frequencies ooj and radii rj are given by 


K-j(r) = -muj]{r -rjf, 


( 3 ) 


with j = i,m,o labeling the inner, middle and outer rings, respectively, as shown in 
Fig. 1. In order to obtain the lowest energy eigenstate of a single atom in a harmonic 
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Figure 1. Sketch of the trapping geometries formed by the combination of concentric 
(a) harmonic and harmonic ring potentials; and (b) two and (c) three harmonic 
ring potentials. Transport will always take place from the innermost trap to the 
outermost one. The top row shows three-dimensional representations of the trapping 
potentials. The bottom row shows 2D schematics of the three configurations, with 
traps represented in gray, and including labels for trap frequencies and radii, ojh is the 
radial trapping frequency of the harmonic potential; rj and uij with j = i,m,o are the 
radii and radial trapping frequencies of the inner, middle, and outer rings, respectively. 


ring potential, we use a variational approach with the ansatz 








( 4 ) 


where aj and fij are the variational parameters, and Nj is fixed by the normalization 
condition 27t = 1, which leads to 




\ 


2a j / TT 


tb + (erf {(djCj^/Zaf) + 1) 


( 5 ) 


The joint potentials VT{r) sketched in Fig. 1 are modeled as piecewise combinations 
of the corresponding potentials given by (1) and (3), and are truncated at the points 
where they coincide (see Eqs. (14), (16) and (17)). Truncated harmonic potentials have 
shown to give qualitatively similar results for SAP as smooth Gaussian [19], Poschl- 
Teller potentials [21], or atom-chip systems [22]. Moreover, they allow to use a simple 
ansatz for the localized ground states such as the ones given in Eqs. (2) and (4). 

The dynamics of a single atom in this potential are then governed by the 2D 
Schrodinger equation 

d 

ih—^{r,(j),t) = H^{r,(j),t), 


where the Hamiltonian is 


( 6 ) 
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Note that the azimuthal dependence (j) in (6) does not play any role in the dynamics if 
they are restricted to the ground states of the harmonic and harmonic ring potentials 
for which 0, f) = f). In the following, we will focus mainly on these states, 

although states carrying 0AM (and thus having an azimuthal phase dependence) will 
be addressed at the end of Section 4. 

Throughout the paper, we will express all variables in dimensionless harmonic 
oscillator (h.o) units with respect to 14 {rn = h = Uh = 1)- 


3. Spatial Adiabatic Passage between two concentric potentials 

3.1. Two-state model 

In this section we introduce the two-state model for the trapping conhgurations shown 
in Figs. 1(a) and (b). Let us assume that at any time the total atomic wavefunction 
can be written as a superposition of the two orthonormalized localized states of each 
individual potential as 

'^2s{r, t) = ai{f)^i{r) ao(t)To(r), (8) 

where a* (oo) is the probability amplitude for the atom to be in the inner (outer) 
potential. 

Introducing the ansatz T 25 into the Schrodinger equation ( 6 ), one obtains the two- 
state model, 

d ( Qi \ ^ ( 0 —J /2 \ f CLi \ 


ih 


dt \ tto 


= h 


-J/2 A 


where 


47r /-oo _ 

J = — / T*(r)Jr\['o(r) rJr, 

n Jo 

‘2'ii f fOO 

A = — ( / \['*(r)JfTo(r) rdr — 

n \Jo 


T*(r)JfTj(r) rdr 


(9) 

(10a) 

(10b) 


are the tunneling rate and energy bias, respectively. The values of J and A used in 
the two-state model simulations are obtained by performing these integrals numerically 
with the corresponding wavefunctions. The states Tj for j = i,o are obtained by 
orthonormalizing the states and Tr,o for the configuration in Fig. 1(a) or and 
o for the configuration in Fig. 1(b). 

The eigenvalues of the standard two-state Hamiltonian in (9) are 


= I (a ± 4A2 + J2) , 


with corresponding eigenstates 

T+(6') = sin(6')Ti — cos(6')To, 
T_(6*) = cos(6')Ti-I- sin(6')To, 
where the mixing angle 9 is given by 

tan(26*) = 


J 

A' 


( 11 ) 

(12a) 

(12b) 

(13) 
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Figure 2. Numerical simulations using the two-state model. Matter-wave RAP 
process to transport a single atom from the localized ground state of a harmonic 
potential to the ground state a harmonic ring one, see Fig. 1(a). As a function of 
time: (a) trapping frequency Wo and (b) radius of the harmonic ring potential Tq; (c) 
tunneling rate J (black solid line) and energy bias A (red dashed line); (d) mixing 
angle 0; (e) eigenvalues E± of the energy eigenstates dt-i. (black solid line) and 
(red dashed line); and (f) populations Pi — |aip (black solid line) and Po = |aop (red 
dashed line) of the harmonic and harmonic ring potentials. All variables are expressed 
in h.o. units with respect to the potential 14,. 


3.2. RAP-like protocol 


By adiabatically following either or a single atom initially loaded into the 
localized ground state of the inner potential can be efficiently and robustly transferred 
to the outer one. Let us assume that initially, at t = to, one has A < 0 with |A| 3> J, 
such that the mixing angle in (13) is given by ^(to) = 7r/2 and, according to (12a), 
^ 25 (to) = To transfer the atom to the outer potential by adiabatically 

following one needs to slowly modify the energy bias and the tunneling rate such 
that at the end of the process, att = t/,A>0 and |A| !:§> J, which, in turn, means 
0(tf) = 0 and '^ 2 s{tf) = T+(t/) = —'^'o- To avoid diabatic transitions and keep the 
atom in the energy eigenstate 'h_|_ for the whole dynamics, it is important to reach 
signihcant tunneling rates when the sign of the energy bias is reversed such that even 
at this moment the energies of the two eigenstates are substantially different, see (11). 

Let us consider the transfer of a single atom trapped in a harmonic potential into 
a harmonic ring, see Fig. 1(a). In this case, the joint potential prohle reads 


VT{r) = 


Vh{r) 

K.o(r) 


for r < rhr, 
for r >fhr, 


(14) 


with Vhr = ToOJoliuJh + uJo). To investigate the SAP process we modify the ring trapping 
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Po{tf) 



200 400 600 800 1000 


Figure 3. Numerical integration of the 2D Schrodinger equation. Final population 
transferred to the harmonic ring trap from the inner harmonic trap as a function of 
the total time of the process and the minimum radius of the ring trap during the 
process. All other parameters are as in Fig. 2. Black lines correspond to contours at 
Po{tf) = 0.99 and the cross indicates the parameter values that were chosen in the 
two-state simulations shown in Fig. 2. All variables are expressed in h.o. units with 
respect to the harmonic potential Vh- 


frequency Uo and radius Tq according to (see Figs. 2(a) and (b)) 

i+(ir 


Uo{t) — CJo + ^0. 




Voit) = Ro + due 2 -^ 


(15a) 


(15b) 


with the parameter values tf = 400, Rq = 4.5, = —1, ooq = 1.7, d^^ = 0.6, a = 8, 

and a = 2/15. These modihcations of the potential yield the temporal variations of the 
energy bias and tunneling shown in Fig. 2(c), calculated with Eqs. (10). 

The temporal evolution of the mixing angle is plotted in Fig. 2(d), while Fig. 2(e) 
shows the evolution of the energy eigenvalues (11), revealing the existence of an energy 
gap during the whole dynamics. The populations of the harmonic and harmonic ring 
potentials are shown in Fig. 2(f), which demonstrate the high fidelity of the transport 
protocol. Let us comment here that, as in all adiabatic protocols, the exact time 
dependence of the parameters is not crucial for the technique to succeed, and even 
withstands some noise [18, 19]. We have chosen this particular time dependence as it 
allows for shorter total times t/ of the protocol with respect to other simpler temporal 
variations of the parameters. 

To check the validity of the two-state model we have also performed numerical 
simulations of the RAP-like protocol using the 2D Schrodinger equation (6) for the 
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Figure 4. Numerical simulations using the two-state model. Matter-wave RAP 
process to transport a single atom between two concentric ring potentials, see Fig. 1(b). 
As a function of time: (a) trapping frequencies uii (black solid line) and Wo (red dashed 
line) of the inner and outer rings; (b) radii (black solid line) and To (red dashed line) 
of the inner and outer rings; (c) tunneling rate J (black solid line) and energy bias 
A (red dashed line); (d) mixing angle 0; (e) eigenvalues E± of the energy eigenstates 
4'+ (black solid line) and 4'_ (red dashed line); (f) populations Pi = |aip (black solid 
line) and Pq = |aop (red dashed line) of the inner and outer rings. All variables are 
expressed in h.o. units with respect to the harmonic potential Vh- 


same parameters as in Fig. 2. Fig. 3 shows the hnal population transferred from the 
harmonic inner potential to the outer ring, Pq, as a function of the total time of the 
process and the minimum radius of the ring during the process. A population transfer 
above 99% is achieved for a broad range of parameter values evidencing the fidelity 
and robustness of the process. The parameters corresponding to the transport protocol 
investigated within the two-state model above are marked with a cross in Fig. 3. 

The previously described RAP-like technique can also be applied to the 
conhguration shown in Fig. 1(b) corresponding to two concentric harmonic ring 
potentials described by the potential 




Vr,i{r) for r < Vio 

Vr,o{r) for r > po 


(16) 


with fio = {viUJi -b roUJo)/{uJi -t- (Uq). As shown in Figs. 4(a) and (b), we change the 
frequency and radius of the outer ring in time following Eqs. (15) with parameters 
Uq = 1.7, 5^ = 0.6, Rq = 7.5, and 5r = 2, while leaving the inner ring fixed with radius 
Vi = 3.5 and frequency Ui = 2. Then, the temporal profiles of the energy bias and the 
tunneling rate, see Fig. 4(c), are such that the mixing angle 9 evolves from tt/2 to 0, see 
Fig. 4(d). In this case, the followed eigenstate 4/_|_ is transformed from initially being the 
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inner ring ground state to finally being the outer ring one. Therefore, the single atom 
is transported from the inner ring to the outer one with very high hdelity, see Fig. 4(f). 


4. Spatial Adiabatic Passage between three concentric potentials 


4 . 1 . Three-state model 

We now consider the trapping conhguration depicted in Fig. 1(c) consisting of three 
concentric ring potentials, constructed from truncated harmonic potentials 

{ Vr,i(r) for r < 

Vr,m(r) for fim<r < frao (17) 

Vr^oir) for r > fmo 

with fim = [ViUJi + r^Um)/{u:i + Um) and fmo = {rmU^m + roUo)/{u}m + UJo)- 

In analogy with the previous section, we assume that the total atomic wavefunction 
can be written at all times as a superposition of the three orthonormalized localized 
ground states of each individual potential as 

'^3s{r, t) = ai{t)^i{r) + a„,{t)^rn{r) + ao{t)^o{r). (18) 


The equations of motion for the probability amplitudes Oj, Om, of the inner, middle, 
and outer ring potentials, respectively, read 



] 


/ 

A, 


0 

\ 

^ ai 

\ 



= h 



0 








1 

0 


Ao 

/ 

\ ao 

/ 


(19) 


where we have taken the energy of the ground state of the middle ring as the energy 
origin, and the dehnitions of the tunnelling rates and energy biases are obtained 
analogously to Eqs. (10) with the corresponding localized wavefunctions. From the 
diagonalization of this Hamiltonian, one obtains three energy eigenstates Tj. and 
with eigenvalues E± and Ed- For Aq = Aj, the so-called spatial dark state, involves 
only the localized ground states of the inner and outer rings. 


^d(0) = cos 0\kj — sin ©To, 


( 20 ) 


with 


tan 0 


Jim! Jmo- 


( 21 ) 


4-2. STIRAP-like protocol 

SAP for a single atom between the innermost and outermost rings can be achieved by 
adiabatically following T^. Starting with the atom in the inner ring, "^ssito) = Trf(to) = 
Tj, and employing a temporal variation of the tunneling rates such that 0, see (21), 
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Figure 5. Numerical simulations using the three-state model. Matter-wave STIRAP 
process to transport a single atom from the innermost to the outermost potential of 
three concentric rings, see Fig. 1(c). As a function of time: (a) radii (black solid 
line), Tm (green dotted line) and Tq (red dashed line) of the inner, middle and outer 
rings; (b) tunneling rates Jim (black solid line) and Jmo (red dashed line); (c) energy 
bias Ai (black solid line) and Aq (red dashed line) of the inner and outer rings with 
respect to the middle one; (d) mixing angle 0; (e) eigenvalues of the energy eigenstates 
(black solid line), '1'^ (green dotted line), and dt- (red dashed line); (f) populations 
Pi = loip (black solid line). Pm = \am\^ (green dotted line), and Po = |aop (red dashed 
line) of the inner, middle, and outer rings. All variables are expressed in h.o. units 
with respect to the harmonic potential Vh- 


slowly changes from 0 to 7r/2, one obtains at the end of the 

process. Note that this temporal variation of the mixing angle means to favor first 
tnnneling between the middle and outer rings and later on, and with an appropriate 
temporal delay, tunneling between the inner and middle rings. 

Within the three-state model, we have investigated the transport of a single atom 
from the localized ground state of the inner ring to the outer one. Fig. 5(a) shows the 
temporal variation of the inner and outer rings radii to achieve the proper variation 
of the tunneling rates, shown in Fig. 5(b). In this case, the time dependences of the 
radii are described as in the previous cases by Eq. (15b), with Rq = 3, 6r = 2 for 
the inner ring and Ro = 12, Sr = —2 for the outer one. The middle ring has a hxed 
radius = 7.5, and, since no detuning is necessary, all ring trapping frequencies are 
kept constant and equal to Ui = ujm = oJq = 2. The parameter setting chosen implies 
that the orthonormal localized ground states are not fully resonant during the whole 
dynamics, see Fig. 5(c). However, since | Aq — Aj| <C Jim, Jmo in the time window where 
tunneling is non-negligible, it is possible to follow from 0 = 0 to 0 = 7r/2 in this 
interval, as depicted in Fig. 5(d). This adiabatic following results in a high hdelity 
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Figure 6. Numerical integration of the 2D Schrodinger equation to simulate the 
matter-wave STIRAP process to transport a single atom from the innermost to the 
outermost potential of three concentric rings, see Fig. 1(c). Time evolution of the 
populations of the inner, middle and outer rings when the initial state (obtained from 
imaginary time evolution) (a) does not posses 0AM or (b) carries one unit of 0AM. 

(c) Snapshots of the 2D atomic probability density for the three times indicated in (a). 

(d) Snapshots of the 2D atomic probability density (top row) and phase distribution 
(bottom row) for the three times indicated in (b). Parameters are the same as in Fig. 5. 


single atom transport from the inner ring to the outer one, with an almost negligible 
population in the middle one for the whole dynamics, see Fig. 5(f). 

Fig. 6(a) shows the temporal evolution of the populations in the inner, middle, and 
outer rings, respectively, obtained from the numerical integration of the corresponding 
2D Schrodinger equation for the same parameter values as used for Fig. 5, and where 
the initial ground state is obtained via imaginary time evolution. The hdelity of the 
SAP process is above 99%. Fig. 6 (c) shows snapshots of the 2D probability density 
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for the three different times indicated in the temporal evolution. The same protocol 
applied to a state with 0AM is plotted in Figs. 6(b) and (d). The initial ground state 
is taken to be the one obtained via imaginary time evolution with an added azimuthal 
phase exp(i0), and thus with a winding number of one. The time evolution of the 
populations from the 2D Schrbdinger equation in Fig. 6(b) shows a final fidelity above 
97%. Moreover, Fig. 6 (d) shows snapshots for the three different times indicated in the 
temporal evolution of both the 2D probability density (top row) and the corresponding 
phase distribution (bottom row). Note, comparing the phase distribution in (d,i) with 
that in (d,iii), that initial and final states have the same 0AM winding number. This 
is because the joint trapping potential consists of three concentric rings and, therefore, 
cylindrical symmetry is preserved during the whole dynamics. 

5. Conclusions 

We have investigated SAP processes for the transport of a single-atom from an inner trap 
to an outer one in three different conhgurations of cylindrically symmetric concentric 
potentials. In all cases, the energy eigenstate followed has been transformed from the 
localized ground state of the inner potential to the localized ground state of the outer 
ring by manipulating the trapping frequencies and/or the radii of the rings. 

Time-averaged adiabatic potentials would be ideal systems to perform these 
techniques because of their strong confinements on the order of kHz, their possibility 
of realising trapping lifetimes on the order of up to one minute [4] and their ability to 
dynamically change the radius of the traps [3] . Our estimations show that the techniques 
presented in this paper can be implemented in these systems with a total process time 
on the order of seconds. Similar total times can be expected for the painted potentials 
|5, 6], 

We have demonstrated the high efficiency of these SAP processes by means of simple 
two- and three-state models. Additionally, we have checked the accuracy of these models 
by comparing their predictions with the numerical integration of the corresponding 2D 
Schrbdinger equation. The latter has also been used to demonstrate that SAP processes 
between cylindrically symmetric concentric potentials can be applied to transport states 
carrying 0AM. 
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